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In this paper we present a theoretical analysis of a recently proposed two-dimensional Cellular 
Automata model for traffic flow in cities with the novel ingredient of turning capability. Numerical 
simulations of this model show that there is a transition between a freely moving phase with high 
velocity to a jammed state with low velocity. We study the dynamics of such a model starting with 
the microscopic evolution equation, which will serve as a basis for further analysis. It is shown 
that a kinetic approach, based on the Boltzmann assumption, is able to provide a reasonably good 
description of the jamming transition. We further introduce a space-time continuous phenomeno- 
logical model leading to a couple of partial differential equations whose preliminary results agree 
rather well with the numerical simulations. 
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Traffic flow in networks is a problem of great technological interest. In systems where good communications are 
crucial (cities, computers, computer networks, etc ...) the understanding of how undesirable, though possible, "traffic" 
jams form and evolve can lead to significant advances in design and planning of optimal strategies [EJ. 
£j[y Traffic flow is a complex collective phenomenon that presents a rich phenomenology and, for the same reason, is 
difficult to deal with, both theoretical and numerically. In spite of that, given its great practical and technological 
impact, traffic car behavior in roads has been extensively studied. Fluid dynamics and kinetic theory have been 
traditionally used || as the theoretical tools to deal with the problem and treatment has been focused mainly on one 
dimensional situations (highways with one or several lanes, with, maybe, some intersections, is the typical example). 
Only recently H models based on cellular automata (CA) have been proposed to study this kind of problems. They 
. present the advantage of being easier to handle numerically while catching the correct behavior of fluid type systems. In 
the last years extensive studies on different one dimensional CA models have been carried out showing good agreement 
with previous theoretical and experimental results. These studies have also revealed new features of traffic flow. It 
has been claimed, for example, the existence of self-organized criticality in the behavior of car traffic in highways [|J. 
The situation is different in two dimensions (2D). The problem has not been treated in the fluid dynamics framework 
and CA models have opened the possibility to study them. 

To our knowledge only two 2D models |]|| have been reported. Both of them present a periodic city with two 
populations of cars moving on it. Ref. [^J is a deterministic model where cars never turn while in the models in Ref. 
[3| are stochastic and cars are allowed to turn with certain probability. Extensive numerical simulations on them j^]- 
\7\ show the presence of a phase transition between a free moving (high velocity) uniform state and a jammed (low 
velocity) state, where the different populations of cars are separated. In the light of recent studies on CA models 
H the arising of this separation can be understood as a consequence of the violation of the semidetailed balance, 
characteristic of traffic models. 

At the moment there is no theoretical basis to explain all these results and it seems neccesary to build one. The 
only exception has been a simple mean- field like study j^] of the jamming trasition of an extension of the model in 
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Ref. §. In this paper we are going to perform a more ambitious theoretical analysis on which further work could be 
based. We will refer, as our working model, to the one called model A in Ref. || (briefly described in ||), but our 
procedure can be easily extended to other 2D CA models. 

We will begin by writing down the microscopic evolution equations for a set of boolean variables that describe 
exactly the dynamics of the system (section 111). These equations should be taken as the starting point for any 
further approximation. Using the microscopic variables we will be able to give an exact expression for the mean 
velocity of the system. We will show in section IV, by means of very simple arguments, how the behavior of the 
mean velocity is correctly predicted in the freely moving phase. In section [V] we will use the tools of Kinetic Theory 
to obtain an approximate macroscopic version of the microscopic equations by introducing an average over a non- 
equilibrium ensemble. We have started by considering the simplest approach: the Boltzmann approximation which 
assumes the breaking of spatial two-point correlations. Though simple this approximation already shows the presence 
of a transition (by predicting the unstability of the uniform phase) and the correct behavior of the velocity before 



and after the transition. It fails however in giving the density at which the transitions occurs. Finally, in section VI 



we present a continuous (in space and time) description of the model. The purpose of this approach is to make the 
connection with the previous knowledge on traffic flow coming from fluid mechanics. The novelty is that, as far as 
we know, this is the first time that a two component fluid model has been presented to study 2D traffic flow. The 
continuous approximation will lead to a system of two coupled nonlinear partial differential equations that govern the 
evolution of the densities of the two types of cars present in the city. In fact they correspond to some limit of the 
Boltzmann equations of section [v[ We will see how these equations also predict a phase transition as well as they 
give the correct behavior of the velocity of the two phases. We will end up by writing our conclusions in section VI] . 



II. MODEL 

As announced in the Introduction, in this paper we are going to deal with model A of Ref. Q . Briefly, this model 
can be described as follows. We have cars moving inside a town. The town is made of one-way perpendicular (L 
horizontal and L vertical) streets arranged in a square lattice with periodic boundary conditions. Vertical streets are 
oriented upwards and horizontal streets are oriented rightwards. There are N cars moving inside this city. Cars sit at 
the crossings, and they can move to one of its nearest neighbors (allowed by the direction of the streets) every time 
step. Two cars cannot be simultaneously at the same crossing. Among the cars, N/2 of them move rightward with 
probability 7 and upward with probability 1 — 7 (for symmetry reasons < 7 < 1/2), and the remaining N/2 do it 
the other way around. Accordingly, half of the cars move mainly upward and half mainly rightward. Finally, there 
are traffic lights that permit vertical motion at even time steps and horizontal motion at odd time steps. 

The dynamics of the model is as follows. Every time step the direction where to move next is decided for every 
car according to its assigned probability. Then, it is checked that the site where it is going to move is empty and 
that the motion is allowed by the traffic light; otherwise the car will not be moved. Finally, all cars that can be 
moved are placed at their destination site and the next time step starts. We want to stress that the whole process is 
carried out simultaneously for all cars, what implies that if a car leaves a site empty at a given time step, t, it will not 
be occupied at the same time step. The fact that traffic lights allow motion alternatively in vertical and horizontal 
streets prevents two cars from colliding at any crossing. 

Notice finally that by setting 7 = we recover the deterministic model of Ref. || . 

The main results of this model are summarized in Fig. which represents a plot of the average velocity (v) vs. the 
density of cars (n) for different values of 7. The most striking feature of this figure is the first-order-like transition 
from a phase in which cars move freely to another in which cars are jammed. We refer to Ref. Q for more details. 
In this reference it was shown that in the jammed phase cars arrange in diagonal strips - with the two types of cars 
roughly separated in two halves, and the question arose whether it was one or more strips. Arguments were given 
supporting the idea that, in very large systems, the jammed phase has got a multistrip structure (in other words, the 
strips of this phase seem to have a "characteristic" size) . In this paper we will try to reproduce analytically the phase 
diagram of Fig. [fl and to provide an answer to what the structure of the jammed phase is like. 



III. MICROSCOPIC EQUATIONS 



The model can be described by a set of boolean variables for every site of the lattice and every time step. This 
variables are the following: 

• nl, the occupation number of site r at time t by a car with probability 7 of moving horizontally (1 if occupied, 
if not); 
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• v\, the occupation number of site r at time t by a car with probability 7 of moving vertically; 

• £*, a random variable which is 1 with probability 7 and with probability 1 — 7; 

• rf T , a random variable which is with probability 7 and 1 with probability 1 — 7; 

• a = t mod 2, which is 1 if traffic lights permit horizontal movements and if they permit vertical movements. 

A given state of the system is completely characterized by specifying the whole set of occupation numbers, with the 
constraint y^ v v\ = (no two cars simultaneously at the same site). The random variables rf r , £*/ and rf T ,) are all 
independent two by two if r 7^ r' or t 7^ i'; furthermore, they are also independent of the occupation numbers and 
of ct*. These random variables are introduced to decide where cars are going to move to: £* = 1 (rf T = 1) decides 
whether a /x(^)-type car placed at r at time t is going to move horizontally and £* = {rf T = 0) vertically. 
With this set of boolean variables, the evolution equations for the occupation numbers turn out to be: 

a4 +1 = a4 { <t'£ + (/4 +x + ^ +x ) 

+ ^+^ r (/4 + y+^ +y ) } (1) 
+ Mr^r { ^fr-xMr-x + 0%-yMr-y }, 

where we have introduced the short-hand notation 6=1 — 6 and where x and y are vectors representing the dis- 
placement of one site to the right and up respectively. The evolution equation for v\ can be obtained from (|lj) by 
simply exchanging fi and v, and £ and r\ (in what follows, whenever we write down an equation for the /x-type cars, 
it must be understood that there is also another one for the ^-type cars, which can be obtained by this or a similar 
replacement). In this way, the complete evolution of the system is specified by the above set of 2L 2 boolean equations 
plus a given initial state {/J%, v®} and periodic boundary conditions. 

Equation (|l]) can be obtained from simple considerations on boolean variables. First of all, each boolean variable 
represents a condition that may or may not hold (for instance, /x' can be read as "if there is a /x-type car at site 
r and time t...", and so on); secondly, the product of boolean variables corresponds to the logical AND operation 
(whose result is always a new boolean variable), and the addition of boolean variables turns out to be the logical OR 
operation (whose result is another boolean variable, provided all terms are mutually exclusive, i.e. of zero product). 
Now the purpose is to build a boolean variable which is 1 if there is a /i-type car at site r at time t+1 (Eq. (|l|)) and 
if there is not. This is simply achieved by considering that (i) if there is a /x-type car at site r and time t there will 
still be there if it cannot move to the randomly chosen neighbouring site (r + x or r + y) , either because it is occupied 
by any car or because the traffic light is red, and (ii) if there is no car at site r and time t there will be one /x-type 
car the next time step if there is one of those cars at any of the neighbouring sites (r — x or r — y) willing to move to 

site r and its traffic light allows the movement. This provides all the terms of Eq. ([!]); for instance, /4c'£ r is 1 only if 
there is a /x-type car at site r at time t, the traffic lights allow horizontal movement but it decides to move vertically; 
Mr^^r^r+x + ^r+x) i s 1 only if there is a /x-type car at site r at time t, the traffic lights allow horizontal movement, 
it decides to move horizontally, but the site r + x is occupied by any car; etc. The meaning of the remaining terms 
can be interpreted in the same way. It can be straightforwarly checked that all the terms are mutually exclusive 
(remember that 66 = 0) and that all possibilities have been considered. 

Now, we can make use of the constraint ix\v * = 0, which lets us write Jl t r T' t r = 1 — /i r — v\, to transform Eq. (|]) into 

A4 +1 = (cr*|* + cr'^K + cr'^*_ x /i*_ x + : *|*_ y M r - y 
+ C*£r/4(Mr+x + ^r+x) ~ C*£r-x/4-x(/4 + l# 

+ O^rMrOr+y + *4 +y ) - O l f r _ y ^ r _ y ((J* + v\) . (2) 

In this form it can be clearly seen that the highest order terms contributing to the evolution equations are quadratic 
in the occupation numbers, something whose relevance will be made clear in the next section. 

The initial configuration must be such that it satisfies the requirements J^r Mr = v r = N/2, and t l r v r — 0, the 
former counts the number of cars of each type whereas the latter express the impossibility of having two cars at the 
same site. It is straightforward to check that 

E^=E^ = f (3) 
r r 

are constants of motion for the evolution equations (Q) and that /x^f* = is maintained in the evolution. 
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We need now a microscopic definition for the average velocity in order to be able to obtain what could be called 
the 'equation of state' (average velocity vs. density of cars) for this model. If we define the mean velocity at time t as 
the number of moved cars divided by ./V, this can be expressed in terms of the occupation numbers: 

^)=2^E{(Mr +1 -M*) 2 +(^ +1 -^) 2 }. (4) 

r 

This equation can be easily visualized if one realizes that — Mr * s zero only if there is no /i-type car at site r at 
times t and t + 1 or if there is one that remains there at both times; in other words, if there is no movement of any 
/j,- type car at site r between t and t+1. Otherwise the value is 1 (if the car arrives at time t + 1) or —1 (if the car 
leaves off at time t + 1). The square removes the sign and the factor 1/2 let the expression count only the cars that 
arrive and not those which leave. The same applies to the term related to i^-type cars. Notice that only one of the 
two terms may be 1. 

The average velocity was computed in Ref. 0] as the limit 

1 T 

{=0 

which, if the system reaches a steady state, corresponds to its mean velocity fj(i|| . 

Before we go any further we would like to remark that all the formalism introduced in this section also applies 
to the deterministic model || by simply setting 7 = 0, although in what follows we will make explicit use of the 
randomness to obtain information out of the evolution equations. 



IV. LOW DENSITY LIMIT 



Figure fi] shows that the behaviour of the v-n curves in the free moving phase is, to a large accuracy, 7-independent 
(for 7 > 0). In fact, they decrease almost linearly with slope —1/2. Reproducing analytically this remarkable property 
is a first important result which can be derived from the microscopic equations. It will be the main goal of this section. 
To this purpose notice that Eq. ^ can be rewritten 



«(t) = i-i([M t+ V] + KV]) , 



(6) 



where the following notation has been introduced: [A 1 ] = L~ 2 ^2 T Al, the average over the whole lattice of a given 
function A l r , and the car density, n, defined as n = NL~ 2 . In getting (g) it has been used that any boolean variable b 
verifies b 2 = b and that [/i 4 ] = [z/] — n/2. Time correlations can be transformed into space correlations by introducing 
the dynamics, i.e., by using the equations of motion (0): 



+ ^{[eV] + [!Vx] + [!V4]} 



(7) 



and the counterpart for the v y s (obtained by exchanging £ and 77 and /1 and v in the equation above). We have denoted 
with the subscripts x or y a translation by those unit vectors. Equation (Q) simplifies in the L — > 00 limit, for then 
the random variables can be replaced by their averages (7 or 1 — 7), as they are independent of the car variables. 
Furthermore, in the low density limit (n — + 0) we can assume that the occupation numbers become independent 
variables in the uniform steady state ([/t*/ix] ~ [/**] [a*x] = n 2 /4, and similarly for the other seven spatial correlations 
- those in (0) and its counterpart. In principle (though we will later come back to this point again) this hypothesis 
is justified by the fact that, due to the random dynamics, at very low densities cars have enough time to "forget" 
between two successive encounters. Therefore, in this regime, 



[M . + vw{,^© 2 } + ^ +2 ,© 2 } 



Plugging this expression into the equation for the mean velocity (g) we get 

w(t)~ i(l-n), 



(8) 



(9) 
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and since this formula corresponds to the uniform steady state it is the average velocity This expression for v is 
7-independent, which is in good agreement with the results obtained in the simulations in the freely moving phase; 
in particular it predicts the slope —1/2 of the v-n curves for small n. 

Clearly, the above argument does not hold for the deterministic case (7 = 0) reported in Ref. M, because in the 
derivation of Eq. (|^) we have made explicit use of the randomness. This is very clear when one compares with the 
result obtained in that reference - v — 1 for the freely moving phase and v — for the jammed phase. However it 
is still an open question whether the jamming transition in this deterministic model drops to n — as the system 
size increases; were this result confirmed Eq. (|9|) would also be valid in the limit 7 — > (remember that its derivation 
assumes that the size of the system is infinite), what would provide a beautifully closed picture of this system. 



V. BOLTZMANN APPROXIMATION 

Equation (|J) and its z^-counterpart make feasible to employ the methods of the standard Kinetic Theory |ll],|l2|] in 
its lattice-gas version ]l3|,0 in order to obtain macroscopic properties of the system. To this purpose let us introduce 
a non-equilibrium ensemble in which all initial conditions, v r }> are equally weighted. We will denote the average 
over this ensemble by (■). When we average over the boolean variables /4 and v\, we obtain real variables - u(r,t) 
and w(r, t) respectively - ranging from to 1 which measure the average number of initial configurations for which 
the site r is occupied at time t by the corresponding type of car (henceforth average occupation). Now, if we average 
over this ensemble the microscopic evolution equations (^]) we obtain 

u(r, t + 1) = (0-7 + CT*7)u(r, t) + a t ^u{r — x,t) + a ju(r — y, t) 
+ o*r/ ((/4/4+x) - (a4-x/4) + (a4^+x) - (/4-x^r}) 

+ CT*7 ((/4/4+y) - <A*r-yMr) + (/4«£fy) ~ (/4-y^)) ■ (10) 

In writing Eq. ( [Tof ) we have used the fact that and r)l are uncorrelated of any combination of occupation numbers 
at the same time step. These equations are not closed due to the presence of the quadratic terms which, on the other 
hand, contain the interaction of the model and therefore cannot be ignored. Evolution equations for the averages of 
these quadratic terms can be written, but cubic and quartic terms will appear, and so on, giving rise to a hierarchy 
similar to the BBGKY one |l2| . Being able to write down a closed set of equations for the average occupations requires 
an approximation. The simplest one is to assume that average occupations at different sites are always uncorrelated. 
In Kinetic Theory this approximation is known as the Boltzmann approximation (or molecular chaos hypothesis). For 
our purpose, this approximation amounts to writing 

(/4/4') = u(r,t)u(r',t) (11a) 
(Mr^') = u{r,t)w(r',t) (lib) 
= w(r,t)w(r',t) (11c) 



for any pair of different sites r and r' and for every time step t. Physically speaking, the molecular chaos hypothesis 
assumes that two colliding "particles" (cars in our model) have never met before, or if they have, all their mutual 
influence has been lost. This hypothesis is expected to be true in the limits of short times (cars do have never met 
before) or vanishing densities (the probability of an encounter is small enough so as to lose all information between 



successive encounters). We have already made use of this fact in Sec. IV. 

A further simplification may be introduced, as it is the replacement of the traffic-light variable, cr*, by its time 
average, 1/2. This approximation is expected to hold in the macroscopic regime (which we are interested in), where 
the microscopic details are not seen. However, the study of microscopic scales would require to keep cr* to its original 
values 0, 1 alternatively. We will go back to this point later. 

In summary, the ensemble average together with these two approximations lead to a closed set of equations (hence- 
forth referred to as the Boltzmann equations) for the averaged variables u(r,t) and w(r,t): 

u(r, t + 1) = -u{t, t) + lu{r -x,t) + |«(r - y, t) 

+ \(r, t){u(r + x, t) + w(r + x, t)} - ^u(r - x, t){u(r, t) + w(r, t)} 

+ |u(r, f){u(r + y, t) + w(r + y, t)} - 2«(r - y, t){u{r, t) + w(r, t)} (12) 
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Iteration of the Boltzmann equations will give a mean-field-like evolution of our model. This evolution will lack some 
of the effects caused by correlations, but it still keeps many of the interesting features of the model. Furthermore, the 
mean field description allows us to obtain analytically many physical qu antit ies of relevant interest. The remaining 



of this section is devoted to study the Boltzmann equations. Subsection V A deals with the linear regime where we 
determine the stability of the homogeneous solution, whereas in subsection |V Bj we present the results of a numerical 
simulation performed directly on the Boltzmann equations and compare them with the linear regime as well as with 
the simulations of the full model Q obtained from the exact microdynamics, Eqs. (||). 

A. Linear Stability Analysis 

The Boltzmann equations ([l2]) have, for any value of the turning parameter 7, the uniform solution 

u(r,t) = w(r,t) = n/2, Vt > . (13) 

This is what simulations show for small enough density However above a threshold (depending of 7) the sim- 
ulations show a steady non-uniform pattern, what means that the uniform solution becomes unstable against small 
perturbations. Thus, settling the stability of ( |l3| ) requires studying the evolution of small perturbations about it, i.e. 
we need to look for solutions of the form 

71 Tl 

u{r,t) = -+6u(v,t), w(r,t) = - + 6w(r,t), (14) 



which introduced into (|12|) give, to linear order: 

$ u (r,t+l) = i{(l + £) Su(r,t) 



Tl 

+ 7— 5u(v + x, t) + 7(1 — n)Su(v — x, t) 

n "1 
+ l^u(r + y,t) +7(1 - n)Su(r-y,t)j 

Tl 

— — {Sw(r 7 1) — jSw(r + x, t) — 7<5w(r + y, t)} ■ (15) 

Equation (^|) is linear but non-local, because it involves the average car occupations at five nodes (r and its four 
nearest neighbours). By taking the discrete Fourier transform, defined as: 

«J/(k,i)=^e- ikr «5/(r,f), (16) 

r 

(where k x ^ — 2Tzq x ^/ 'L, q x (y) = 0, 1, 2, . . . , L — 1 and / is either fi or v) Eq. ( j^ ) becomes 

(Sr + io^( k >(«(M))- <i7) 

where fJ(k) is the k-dependent 2x2 linear evolution operator whose elements are given by 

I 77 1 — 77 

fin = 2 + 4 i 1 + IS* + fS y } + ^-{7^ + IS;} , (18a) 

Tl 

n 12 = -{ 7 s x + is y -i}, (i8b) 

with S x ( y ) = e 4kx («' . The elements O21 and Q22 are obtained from Eqs. ([l^) by exchanging 7 and 7 in SI12 and fin 
respectively. Iteration of Eq. ( |T7| ) yields the dynamics of the system at the linear level in Boltzmann approximation. 
The information on the stability of the uniform state is contained in the eigenvalues of the linear evolution operator 
f2(k), written for convenience as 

n(k)^-(k) = e»'< k >Vi(k), 3 = 1,2. (19) 
Hence the dynamics at the linear level is simply expressed as 
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(LT4)^(<Ce,«0^< k K^ k 4)- (20) 

where the first (second) column of U(k) is the eigenvector V"i(2) s -) ■ Note here that the two eigenvalues must satisfy 

limz,(k)=0 (21) 

as a consequence of the conservation laws (^|) (for the sake of simplicity we will be referring to Zj(k) as the eigevalues, 
though they are actually their logarithm). 

Whenever Re Zj(k) < (j = 1,2) for all k, the amplitude of all eigenmodes, V'jO*), will vanish exponentially 
with time and finally disappear, recovering the uniform state. However, as soon as a value of k appears such that 
Re Zj m (k) > (for j m = 1 or 2, or both), then the amplitude of the corresponding eigenmode ipj m (k) will grow 
exponentially with time, in other words, the uniform state will become unstable, and inhomogeneous spatial structure 
will develop. Therefore, the spectrum of the linear evolution operator determines the macroscopic behavior of the 
system. There is more information that can be extracted from the spectrum when the uniform state is unstable [ p"5| . 
In general, there will be a wavevector, k m , for which Re z Jm (k m ) (for some j m = 1,2) is maximum. Fluctuations 
with this wavevector will dominate as t — > oo and consequently, the structure will possess a typical lengthscale of 
A rn = 2tt /k m , and will be formed along the direction of k m . Furthermore, the onset time of the instability is given by 



[Rez Jm (k m )}- L (see Eq. (|20|)). 

The eigenvalues Zj(k) are easily obtained from the elements of the linear evolution operator, given in (|l8|). A typical 
plot of Re Zj m (k) is shown in Fig. || for n = 0.8 and 7 = 0.1. According to the numerical simulations (Fig. H) for 
this parameter set the system is jammed, i.e., the system is non-uniform since cars cluster in bands which form an 
angle 8 — 7r/4 with the positive x-axis. From Fig. g we see that the linear stability analysis agrees with this fact as 
there exists regions where Re z Jm (k) is positive with maxima located on a line forming an angle of 9 = 37r/4 with the 
positive a;— axis, i.e., perpendicular to the bands of the system. From the location of the maxima we can obtain the 
typical size of the bands or, equivalently, the number of bands which would appear in a finite size simulation. We will 
give more details about this point later in this section. We have not shown the second eigenvalue, because it satisfies 
Re z(k) < for all values of k, so it does not affect the previous discussion. 

The traffic jams of this model are always bands as those of Fig. [j]b (i.e. with 6 = 7r/4). Hence, the maxima of 
Re Zi m (k) appear along the direction given by k = (—1/ \/2, 1 / \/2) . Therefore, in order to look for these maxima in the 
expressions for Zi(k) is enough to choose a wavevector of the form k = /ck. This means that S x = S* = exp(— ifc/v2). 
For this particular direction the eigenvalues take a simple form: 

e*±< fc > =1-Wl-- 



2 

± s |£^! - 2(1 - n) Q - n) C 2 (l - s 2 )} 1 ' , (22) 

where s = 8^(^/2-^/2) and ( = 1 — 27. It is straightforward to check that Re z+{k) > 1 for some /c's if and only if 
n > 1/2. In consequence, this analysis predicts a jamming transition for n > 1/2 irrespective the value of 7. The 
existence of this transition is the most important result of this section; furthermore, the prediction for the transition 
density is reasonably good for 7 close to 1/2, the highest randomness (see Fig. |l|a). However it disagrees with the 
simulations for smaller 7's. 

A comment about the role of the traffic light variables, cr* is needed here. Figure || (and the discussion about it) 
has been obtained for er* = 1/2. Keeping the variables er* to their original values, 0, 1, alternatively, one obtains two 
evolution operators, one for even time steps, f2 e (k), and another one for odd time steps, tt° (k). Then the eigenvalue 
problem in Eq. (|lj]) has to be formulated in terms of an effective evolution operator given by {f2 e (k)ri°(k)} 1 / 2 . 
However, in doing so, one has to face the fact that f2°(k) and f2 e (k) do not commute. Thus we end up with three 
sets of eigenvalues, those of Eq. (|lj|, of {r2 e (k)fi°(k)} 1 / 2 and of {fi°(k)fi e (k)} 1 / 2 . The largest difference between 
these three sets is for the wavevectors k located at the edges of the Brillouin zone (k ~ tt) , related to the microscopic 
behaviour, when the system is observed at scales of the order of a few lattice spacings. On the contrary, for wavevectors 
at the center of the Brillouin zone (k ~ 0) the three sets coincide. This is the region we are interested in, because 
it describes the macroscopic behaviour. As a summary, the macroscopic behavior is not affected by the traffic lights 
being 0,1 or its average 1/2. For this reason, and for the sake of simplicity, we have made all the calculations with er* 
replaced by 1/2. 

According to Kinetic Theory, for each quantity conserved by the dynamics there exists an eigenvector as given by 
Eq. (|l^) whose eigenvalue satisfies the relation (|2l|). In the long wavelength (small wavenumber) limit, the eigenvalues 
can be written as: 
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Zj-(k) = ikc j + (ik) 2 D J + --- (k^O), (23) 

The quantity Cj is interpreted as the speed of propagation of the mode, while Dj is the corresponding diffusivity. In 
the case of fluid dynamics in three dimensions, there are five conserved quantities, and five of such modes: two sound 
modes, two shear modes, and finally a thermal mode. This is a general result of Kinetic Theory and also holds for 
lattice gases as models for thermal fluids In the case of a model of concentration diffusion leading to the Fick 
law, there is only one mode (the concentration); cj is then the drift velocity and Dj the diffussion coefficient. The two 
eigenvalues of our model (E g. (p2|) ) correspond to the two conservation laws (|^), so that expanding them in powers of 
the wavevector k, as in Eq. (|23|), we can identify the two (opposite) speeds of propagation and the (unique) diffusion 
coefficient, with the result (for n < 1/2): 

1 / 1 \ 1/2 

c± = ±-((l-n) l / 2 (n-~) (24a) 

D = [2 - n - (1 - n) (1 - 2n)C 2 ] (24b) 

We have to state here that this velocity has nothing to do with the remnant velocity measured from the simulations, but 
it is the velocity with which long-wavelenght fluctuations propagate (the so-called kinematic velocity in fluid-dynamics 
studies of traffic flow S). 

As it was said before, by calculating the maximum with respect to k of Re z+(k), with k along the direction of 
9 = Sit/A (Eq. (p2f)) we obtain the characteristic length of the system as A m = 2ir jk m . The relevance of this quantity is 
that it provides information about the typical center-to-center band separation in the jammed phase and consequently 
determines the number of bands. We will not write here the expression for A m , because it is very cumbersome. Instead, 
in Fig. |^ we plot the value of A m as a function of the density n for some values of the randomness 7. It can be seen 
that apart from both edges (n > 1/2 and n < 1), the bands have a typical width of few tens of lattice spacings, 
in qualitative agreement with the simulations. However, in actual simulations the characteristic length may not be 
the natural one (that the system would have if it were infinite) because of a competition with the finite size of the 
system. If the system is not large enough to form an integer number of bands of width A m it will form less bands, but 
concentrating all the cars in the existing bands. Therefore the characteristic length will be modified. However, the 
value of X m we have obtained can still give a good prediction of the most stable number of bands in the simulations 
at finite L. The number of bands predicted by the Boltzmann theory is simply L/X m . If L/X m is not an integer, the 
system will have the trend to form \_L/X m \ bands (where \_x\ means the largest integer smaller or equal to x). Reading 
off the Fig. H for the case 7 = 0.4 we obtain \L/X m \ — 1 for all n > 1/2, in perfect agreement with what is shown Fig. 
|l|a (made for a city with L = 64). The absence ofjumps in this curve (in the jammed phase) means that only one 
band has been formed. For the case 7 = 0.3, Fig. | predicts [L/X m \ = 1 for n < 0.6, [L/X m \ = 2 for 0.6 < n < 0.7 
and [L/X m \ = 3 for 0.7 < n < 0.95. Inspection of Fig. [I] shows that there is only one band for n < 0.65, while for 
n > 0.65 there are two bands, as indicated from the increase of the velocity. For higher densities, although it cannot 
be seen in Fig. |l], occasionally some simulations end up in a three- band structure ||. Unfortunaltely, the agreement at 
this level between the Boltzmann predictions and the simulations is poorer the smaller the randomness. For 7 = 0.2 
and 7 = 0.1 the number of bands given by the Boltzmann theory for a system size of L — 64 at intermediate densities 
is about 5 or 6, while in actual simulations such number is smaller (1, 2 or 3 bands, but never more). This result could 
be explained by the fact that the higher the number of bands the more difficult to see them in simulations, because of 
band coalescence at their early stages of formation. On the other hand, the Boltzmann approximation was shown to 
be worse the lower 7, because it predicted the jamming transition for n = 1/2, a result that can be valid for 7 close 
to 1/2 (see Fig. [I]), but certainly is not for small 7, where experimentally one finds that it can be as low as n ~ 0.24. 

At both edges of the plot presented in Fig. || the quantity A m goes to infinity as A m ~ [n — 1/2) -1 / 2 for n > 0.5 
while A m ~ (1 — rt) -1 / 4 when n < 1. The former suggests that just after the transition only one band is formed; the 
latter means that all bands coalesce into one which fills up the whole system. As X m — > 00, k m — > 0; thus, as the 
eigenvalues obey the restriction (^l|), the onset time of the instability (given by [Re z(k m )] _1 ) diverges just above the 
transition. This is indeed observed in the numerical simulations, where extremely long relaxation times are needed 
near the jamming transition. For 7's close to 1/2 the onset time of the instability also diverges (see Eq. (p2[)). For 
intermediate n (in the jammed phase) the typical onset times range from a few hundreds to a few thousands time 
steps. 

The last remark of this subsection concerns the special case 7 = 1/2. For this value, for which there is no distinction 
between cars of type fj, and u, Eq. ( p2] ) predicts negative eigenvalues for any value of the density. As a consequence, 
there is no jamming transition for this particular case, in agreement with the numerical simulations (Fig. ^a). 
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B. Simulation of the Boltzmann Equations 



All the information presented in the previous subsection was obtained within the Boltzmann approximation, but 
in the linear regime. This regime is only valid for small perturbations around the homogeneous state. However when 
the perturbations are not longer small, as it happens once the jammed phase appears, the full nonlinear Boltzmann 
equation takes over, and the information given by the linear theory may not hold. In particular, no linear theory can 
give the saturation (i.e. completely filled nodes) occurring in the bands in the jammed phase. 

In order to study the jammed phase we have performed simulations of the full Boltzmann equations (|l^), starting 
with a uniform state slightly perturbed at random. The evolution equations ( |l2|) are then iterated until a stationary 
state is reached. This final state can be either again uniform (when it is stable) or non-uniform, exhibiting a traffic 
jam. This latter situation is illustrated by Fig. |f| for n = 0.6 and 7 = 0.2, where a plot of the average occupation 
profiles in the steady state is shown on the line perpendicular to the bands. We clearly see the band structure, with 
a completely saturated region (u = 1 or w = 1), bounded by partly filled, very narrow layers. The result is similar 
to that obtained in direct simulations on the microscopic model ||. The center-to-center band distance, as measured 
from the simulations is equal to 17 (in lattice units), while the linear stability analysis predicts A m ~ 16, in very good 
agreement. 

The phase transition is numerically obtained for the same region of parameters as those predicted by the linear 
theory (n — 1/2 for all 7). This result supports the idea that the fact that the transition is 7-independent is an effect 
of the Boltzmann approximation itself, and not of the linearization. In fact, from the general theory of stability of 
dymanical systems, one can prove that the stability is given by the linearized equations (unless zero eigenvalues are 
obtained), not being modified by the nonlinearities. Therefore it is the breakup of the correlations what gives the 
degeneracy in the transition density. So, finding a 7-dependent transition requires going beyond Boltzmann. Some 
work along this line is currently being done. 

From the simulation of the full Boltzmann equations we can obtain a quantitative prediction for the velocity in the 
jammed phase. In order to measure the velocity we cannot use the microscopic definition given in Eq. (^), because 
in the steady state the average occupations do not change in time. However, we can use the modified Eq. (R), where 
no time variations appear, to compute the velocity. Figure [5] summarizes these results. In spite of the 7-independent 
transition, there is a strong similarity between this figure and the phase diagram obtained from the simulations of 
the microscopic model (Fig. Q); in particular, the slope —1/2 before the transition is properly obtained, as well as the 
behavior of the remnant velocities after the transition has taken place. 

In summary the Boltzmann approximation gives good predictions for many physical quantities of the model, as 
remnant velocity, size of the bands, onset time of instabilities, and predicts the existence of a jamming transition, 
although it does not capture the fact that the transition is 7-dependent. But with this simple scheme one can obtain 
a qualitatively correct picture of the behaviour of the system. 



VI. A PHENOMENOLOGICAL MODEL 

We have seen in Sees. Ill and 0, respectively, an exact microscopic and an approximate Boltzmann descriptions 
of the model. Both of them consider space-time variables that are discrete. In this section we are going to present 
a continuous phenomelogical approach. Our motivation is twofold. In the first place it provides a new theoretical 
tool to understand the dynamics of the system. In the second place it makes a connection with the most usual 
approaches to traffic flow problems: continuous fluid dynamics [Q. This section only tries to fix the main lines of a 
more complete and exhaustive study of this traffic model from the point of view of continuous dynamics that will be 
presented elsewhere. 



A. Equations 



To define the continuous model we have to consider first the important scales in the system, then redefine the 
variables, and finally take the appropiate limit. Let us assume that the distance between city crossings is e, that the 
linear size of the system is C = Le and that there is a scale of speed for the cars c, defined as the speed of a car that 
moves freely (that is, the length it will travel, if it does not find any red light or other cars, divided by the spent 
time). The time interval taken by half a cycle of a traffic light (one time step) is then e/c. 

Let us define a new set of space-time variables: x = n x 5, y = n y S and r = n T 6/c, where S = le <C 1 is a coarse 
grain length scale (I is an integer such that 1 <C I <C L) we will define shortly, and n Xl n y , n T are positive integers in 







the interval [0, L/l}. Now we introduce, as the variables of our model, two density functions, one for each type of car: 
u(x,y;r) and w(x,y;r) 1 representing block averages in space and time of the boolean variables ^* and v\: 

u(x,y;r) = ^ ^ a4> (25) 

(r,f)eC (l)5) x7 T 

w(x,y 7 T) = ^ Y ( 26 ) 

(r,f)€C (l)B) x7 T 

where Cr XiV \ is a square of size 6x5 (I ) sites, centered at (x, y) = (n x S, n y S) and I T is a time interval consisting of I 
traffic light half-cycles (i.e. it lasts a time St = 5/c). Notice that by definition the densities u and w are positive and 
less than one: 

0<u{x,y,T) , w(x,y,T) < 1, (27) 

and the dynamics, through the impossibility of having two cars at the same site, imposses: 

S(x,y,T) < 1, (28) 

where we have defined the new variable S(x, y; r) = u{x 1 y; t) + w(x, y; t), that stands for the total density of cars. 

Now let us take the continuous limit by making n x ,n y ,n T — > oo, e — > and l,L — > oo while keeping <5 — > and 
(x,y) gfi, r e 7\L + , and £ finite, where: 

ft = {(z,y) : < a;,y < £} e "ft 2 . 

In this way the functions u(x,y;r) and w(x,y;r) represent coarse-grained local densities of the two types of cars. 
Because of this definition these density functions lose track of the microscopic details of the model (at scales e), and 
they are supposed to be almost constant up to scales of order S. 

Let us now write the equations that govern the evolution of u(x, y; r) and w(x, y. r). We do so by writing a flux- 
balance equation for the cars of each type going in and out of a block Cr x ^ y y that contains I 2 ^> 1 sites, with densities 
(constant by definition) u(x, y; r) and w(x, y; t). The cars going in are located either in the block Ct x _g tV y inmcdiately 
to the left, or in the block C( x . y _$) inmediately below. The cars going out will go to the blocks C( x+ g iV y located to 
the right, or to the one above, C( XtV+ gy Let us estimate how many cars enter C( x . y ) from C^ x s, y ) in a time interval 
S t . In the first traffic light half-cycle only cars in the column sited at the the right border of Ci x _g iy \ have any chance 
to go into Ct XtV y In average, the number of cars of both types on that column are u(x — S, y; r)l and w(x — 5, y. r)l. 
These cars have the chance to cross the left border of C^y) and sit in its leftmost column. This column is occupied, in 
average, by S(x, y; t)1 cars, leaving [1 — S(x, y; t)]1 empty sites. The number of cars that actually go through depends 
on the relative positions of the cars in both sides of the line and on the randomness built into the car movements. For 
simplicity we assume that the distribution of cars on both sides of the border is the most probable one: that in which 
the cars are uniformly distributed along the line of length 6. In this case we can assume that in the first time step of 
St, ju(x — S, y. t)[1 — S(x, y; t)]1 cars of type /i and ^w(x — S, y; r)[l — S(x, y; t)]1 of type v will go from C( x -s.y) to 
Cr XtV y During this first time step there has been also movements of cars inside C/ x _g iV y but we are assuming that 
the densities u(x,y;r) and w(x,y;r) are constant in space, throughout the whole element, and in a time interval St- 
We will see later that this assumption is consistent with the resulting equations. Under this condition we can suppose 
that the rightmost column of Ct x _$^ y y and the leftmost column of Cr x<y y are always occupied by the same number of 
cars. Therefore the total number of cars of type u that will cross during S t (I time steps) the left border of C( Xi!/ ) is: 

^ = 7 u(x - S, y; r) I [1 - S(x, y;r)}l. (29) 

In the same way, the number of cars of type u that will go inside C( Xty ) from the block below, C( x>y _$y is 

N? n T = 7 u(x, v -5;t)1[1- S(x, y r)] I, (30) 

with the accompanying equations for iV™ _^ and ./V™ ^ exchanging u and w and 7 and 7. Similarly we can obtain 
expressions for the number of cars going out of C^ x y y Putting all together we can write the balance equation for the 
increment in the number of cars of each type in Cf X)V y. 

2l 2 [u(x, y;r + S t ) - u(x, y; r)] = 

— in — * in ] > out t ou ^ 

= I 2 [711(1 - S,y;r) +ju(x,y - S;t)][1 - S(x,y;r)] - 

- I 2 u(x, y t) [1 - 7 S{x + S, y r) - 7 S(x, y + S; r)] (31) 
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and the corresponding expression for w. The factor 2 comes from considering that only cars from one of the sides 
can enter Ci x ^ v \ at each traffic light half-cycle. Observe that if we suppose that the spatial variables in Eq. (|l^) are 
continuous, we get exactly Eqs. (|3l|). In this way we can think of the coarse-grain treatment in this section as a 
continuous version of the Bolztmann approximation of section Taylor expanding Eqs. ([^) to first order in 5 and 
5t and fixing c = 1 to be consistent with the discrete model, we get the equation for the time evolution of u(x,y;r) 
and w(x, y; r): 

Ur = \{l [u(S - 1)} X + ^{u(S - l)] y } (32a) 

wr = ^{l[w(S - 1)} X + j[w(S - 1)U (32b) 

where the subscripts r, x and y mean, respectively, d/dr, d/dx and d/dy. These equations have to be completed with 
the periodic boundary conditions in C and initial conditions. 

If Eqs. (|52|) represent faithfully the original microscopic discrete traffic flow model they should fulfill a few neccesary 
(though, of course, not sufficient) conditions: conservation of the number of cars, positivy of the densities, excluded 
volume and the existence of stationary solutions describing the two observed phases in the simulations. Let us check 
all of them. First observe that Eqs. (|2|) satisfy the following conservation laws (constant number of cars of each type) 
for any r € 1Z + : 



u(x,y;r)dxdy — / w(x,y,T)dxdy = nC /2 . (33) 
si Jn 

Second, both densities of cars have to stay positive during the evolution if the initial conditions were positive, and 
also the total number of cars per site should be less than one (excluded volume) (^8|). It can be shown from a 
maximum principle |]l9f that both restrictions hold during the evolution if the initial data u(x,y;0) and w(x,y;0) 
satisfy them. And finally, from numerical simulations it is expected that the density functions u and w corresponding 
to a stationary uniform state will be solutions of equations (B2) . It is also expected that a function of band type (see 
Fig. [l]b), representing the jammed state, will be a solution. It can be inmediately seen that the uniform solution: 

u(x,y,r) = n/2, 

w(x,y;T) = n/2, (34) 

is stationary. It can also be checked that a jammed band configuration (of course parallel multiple bands are also 
possible): 

u(x, y; t ) = 9(x - x - y)9(a + y - x + x ), 

w(x, y; t ) = 9(x - x + a - y)8(y - x + x ), (35) 

where 9(x) is the Heaviside step function, is a stationary solution too, but in the sense of the distributions ppj . 
Observe that the arguments of 9 should be taken mod C due to the periodicity of the city. The arbitrary parameter 
xq fixes the position of the band and 2a = nC is its width. 



B. Stability 



We will proceed now to perform a stability analysis similar to that in section V A . To study the linear stability of 
the solution ( |34| ) let us perturb around the uniform solution: 

u(x, y; t) = n/2 + Su(x, y; r), (36a) 
w(x,y;r) = n/2 + Sw(x,y;r), (36b) 

where Su, Sw -C n. Substituing in ( |3^ ) we get, to first order: 

2P T = AP X + BP y , (37) 



where 



Sw) \ 7f 7(§»-l) 



(38) 
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and the expression for B is the same as the one for A exchanging 7 and 7. Expanding in its Fourier series, Eq. (|3 
can be easily integrated. If we define the Fourier coefficients: 

P{k x ,k y ;r) = f ^ P ( x+ , x _ ;r )e- ik * x - ik «y, (39) 
Jn L - 

with fc xfe) = 2Tm x{y) /C,n x{y) £ Z, (H) becomes: 

2P T = ik x AP + ik v BP, (40) 

whose solution is: 

P(k x ,k y ;r) = e 4 ^ A+fc * B > P(k x ,k y ;0). (41) 

The eigenvalues of k x A + k y B are: 



Ax, 2 = k + -n-l)± J 2(^(1 - n ) - - » + " (42) 



where k± = (k x ± k y )/2. The homogenous solution (|3J) will get unstable if the term inside the square root in (|4 
becomes negative, that is, if 

2C 2 ^(l-n)Q-n)+fc^<0 (43) 

If n > g condition ([l3|) is always fulfilled, for any value of k + , for some big enough value of fc_. Once again we obtain 
that the uniform solution becomes unstable for n > i. Observe that Eq. (43) agrees to first order with the result 
given in Eq. ( |22j ) when we make (k x , k y ) = (— fc/\/2, kj\f 7 £). This is because Eq. ( |4C| ) is just the continuous version of 
Eq. (|l7]) to order k, showing again the correspondence between both descriptions(see also the comments below Eq. 
©)■ 



C. Velocity 



Following the reasoning of Sec. VIA we can also get a phenomcnological estimate of the instantaneous average 
velocity in this model. 

The velocity is proportional to the number of movements that take place in the system at any given time. We can 
measure the total number of produced movements by estimating how many cars cross horizontal and vertical line 
elements of lenght S. For example the number of cars of type u and w crossing a vertical line of length 5 centered 
at (x,y) in a time step is (see Eq. ( p9| ) for a similar reasoning), respectively, u(x — e/2,y;r)/[l — S(x + e/2,y;r)]7 
and w(x — e/2,?/;r)Z[l — S(x + e/2,y;r)]7 where we have explicitly pointed out that the cars crossing the line are 
just half a site to its left and they move to a place half a site to its right, and therefore the corresponding densities 
should be evaluated at x — e/2 and x + e/2, respectively. Similarly, the number of cars of type u and w crossing 
a horizontal line of length 6 centered at (x,y) is (see (3^)), respectively, u(x,y — e/2;r)l[l — S(x,y + e/2;r)]7 and 
w(x, y — e/2; t)1[1 — S(x, y + e/2; t)]^. But when we take the limit e — > we can replace u(x ± e/2, y\ r) by 

u{x ± 1 y- ) r)= lim u(x ± e/2,y;r), (44) 

and u(x, y ± e/2; r) by 

u(x,y ± ;T)= lim u(x, y ± e/2;r). (45) 

e-»0+ 

Taking care of all the terms in the same way, properly normalizing (the factor 1/2 comes from the traffic lights) and 
taking the limit e — > 0, L — > 00, C finite, we get the total speed of the system as: 



1 



dxdy{[\ - S(x,y + ;T)][ju(x,y ;T)+^w(x,y ;r)] 



n 



2nC 2 

+ [1 - S{x + ,y;T)}[ju(x^,y,T)} +jw(x^,y;r)}. (46) 
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With this expression we can easily calculate the speed of the homogenous solution ( |34| ) 

1 — n , 
v = (47) 

in accord with all the previous models and the simulations. We can also calculate the speed for a solution of the 
jammed band type ( |35| ) getting, before taking the limit inside the integral: 

V = -tt-c = lim (48) 
Cn L^oo Ln 

When we take the limit e — > (keeping £ and n finite) we get v = 0. This result can be easily understood. In one 
band the only possible movements are along its border (~ jL) and therefore the contribution to the velocity is of 
order ^jL/nLi 2 . This is exactly Eq. ([48j). This makes us to conclude that an infinite number of bands is neccesary to 
have nonzero velocity in the jammed phase of an infinite system. For example, if there is a band every A sites (see 
the discussion below Eq. (pl|)), the velocity will be of order 

JJl = JL (49) 
nL A n y [ ' 

which is finite in the L — > oo limit. 

Let us return to Eq. ([48]). Though it is zero in the e — ► limit it gives a prediction on the behavior of the v — n 
curves after the transition when L is finite and e > 0, 

V = h ^ 

To compare this prediction with the results of simulations we have taken data from Fig. [|. We have considered only 
the points that represent stationary final states with one band 6[ and we have plotted in Fig. ^| the ratio R = Lnv/j 
against n for four different values of 7. If the prediction of Eq. (50|) were completely correct the four graphs should be 
the horizontal lines R = 1. Instead, as can be seen in the figure, the data fit well to four horizontal lines with different 
values of R. The first important conclusion to draw from Fig. ^ is that, because the data fit to horizontal lines, Eq. 
( pp| ) gives the correct behavior of the velocity after the transition: decreases with n as 1/n and is proportional to 
7/L. The discrepancy is in the proportionality factor, and may be due to the fact that the jammed phase obtained 
in the simulations is not a "discontinuous" band (that goes sharply from one to zero in a distance e) as the one in 
Eq. (|3^). Observing carefully the final configurations it can be seen that in fact the intermeditate region where the 
density goes from one to zero, though narrow, extends for several lattice spacings, being narrower the smaller is 7 
(the system is less random). The values of R in Fig. ^ confirm this behavior too. 

VII. CONCLUSIONS 

In this paper we have developed a theoretical framework to describe the two-dimensional cellular automata traffic 
flow model introduced in Ref. The starting point is the exact microscopic evolution equations (|l|) together with 
the microscopic definition of the mean velocity (Q). From these basic ingredients several approximate schemes can 
be carried out. The simplest one is to follow the standard kinetic recipe of breaking up spatial correlations into 
product of one particle averages. This mean-field-like approach leads to a reasonably good description of the model. 
In the first place, the 7-independent (for 7 > 0) slope of the v — n curves in the freely moving phase is correctly 
obtained. In the second place, the existence of a phase transition between the freely moving and the jammed states 
is also explained by such a simple approximation. Nevertheless, it fails to describe properly this transition because 
it predicts a 7-independent transition density. Finally, this theory allows us to understand the rich structure of the 
jammed phase. In this phase, the Boltzmann approximation is able to explain the formation of bands completely 
filled with cars, and the phenomenon! of phase separation of the two populations of cars inside these bands. Besides, 
it gives an answer to the question posed by the simulations about the number of bands in the jammed state. Our 
Boltzmann theory gives an estimation of the average band-to-band distance. This distance appears to be finite for all 
densities in the range 1/2 < n < 1 (excluding the transition density according to this approximation, n = 1/2, and 
the full system n = 1). The conclusion is that there is an infinite number of bands in an infinite system. On the other 
hand, the band-to-band distance estimation let us evaluate the most probable number of bands for finite L, once the 
values of 7 and n are fixed. The obtained results are in good agreement with the simulations. 



In Sec. VI we have faced the description of the model with an utterly different phenomenological approach. This 



approach ends up in a continuous model whose evolution is governed by a couple partial differential equations. 
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Although this approach seems to have no relation to the Boltzmann approach, we have shown that there is indeed a 
connection to the kinetic formalism - from which it can be obtained in the large spatial scale regime. The continuous 
approximation leads to a connection between our model and those of fluid dynamics - the most usual tools to study 
one-dimensional traffic-flow problems. Apart from its theoretical relevance, this phenomenological model is able to 
provide a 7-dependent expression for the velocity as a function of n in the jammed phase. This result, together with 
the linearly decreasing velocity obtained for the freely moving phase, completes a theoretical description of the v-n 
phase diagram. 

The theory presented in this paper for a particular model opens new ways in the theoretical study of two-dimensional 
traffic flow, which deserve future development. First of all, it should be clear that equations similar to (|l|) can be 
written for any cellular-automata-based traffic model and, therefore, the same approximations can be applied to get 
information out of it. ft is also feasible to write down a continuous model, either by itself or as a limit of the kinetic 
description. Secondly, it is possible to devise another approximations based on the microscopic equations ([l]). For 
instance, it is our belief that the next step in this scheme - namely including two-point spatial correlations - might 
account for the 7 dependence of the transition density. We are presently working on it. Finally, it is worth to remark 
that the approximate models reported in this work are far from being exhausted. Questions such as the mathematical 
properties of the instabilities in the continuous model or how to get an evolution equation for the velocity remain 
unsolved. 
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FIG. 1. Results of the simulations for the model of Ref. || for a city with L — 64: (a) phase diagram of the model - average 
velocity (v) vs. car density (n) -, exhibiting the transitions from the freely moving phase to the jammed phase for a set of 
values of 7 (see Q] for details); (b) snapshot of a final run with n — 0.672 and 7 = 0.2 (white and grey represent the two types 
of cars, while black represents emty sites). Notice the phase separation of the two populations of cars and the band- like traffic 
jam. 

FIG. 2. Plot of the real part of the largest eigenvalue of the linear Boltzmann operator (only positive parts are seen). The set 
of parameters is n — 0.8 and 7 = 0.1. Regions where Re z(k) is positive appear along the direction forming an angle 6 = 37r/4 
with the a>axis (see text). 



FIG. 3. Typical band-to-band distance, A m , in the jammed phase, vs. the density, n, for several values of the randomness 7. 

FIG. 4. Jammed phase for the set of parameters n — 0.6, 7 = 0.2, as obtained from the simulation of the full Boltzmann 
equations (jl^) on a lattice of size 64x64. We only plot here the occupation profiles along the line perpendicular to the band. 
Cars of type fi (y) are represented by a dotted (dashed) line. The total population if plotted with a solid line. 



FIG. 5. Average velocity of cars, v, as obtained from the simulation of the non-linear Boltzmann equations (|12]), vs. the 
car density, n. The transition to the jammed phase occurs at n = 1/2, regardless of the value of the turning parameter 7, in 
contrast with the microscopic simulations (compare with Fig. hi). 



FIG. 6. R = vLn/'y against the car density, n, for the data taken from Fig. hj in the single-band jammed phase (see discussion 
after Eq. 
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